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Abstract 



We study the optimality of linear precoding for the two-receiver multiple-input multiple-output 

(MIMO) Gaussian broadcast channel (BC) with confidential messages. Secret dirty-paper coding (S- 

\Q DPC) is optimal under an input covariance constraint, but there is no computable secrecy capacity 

y—( expression for the general MIMO case under an average power constraint. In principle, for this case, 

the secrecy capacity region could be found through an exhaustive search over the set of all possible 

r_^ matrix power constraints. Clearly, this search, coupled with the complexity of dirty-paper encoding 

I— | and decoding, motivates the consideration of low complexity linear precoding as an alternative. We 

prove that for a two-user MIMO Gaussian BC under an input covariance constraint, linear precoding is 

optimal and achieves the same secrecy rate region as S-DPC if the input covariance constraint satisfies 

a specific condition, and we characterize the corresponding optimal linear precoders. We then use this 

t— I result to derive a closed-form sub-optimal algorithm based on linear precoding for an average power 

constraint. Numerical results indicate that the secrecy rate region achieved by this algorithm is close to 

that obtained by the optimal S-DPC approach with a search over all suitable input covariance matrices. 

^£ I. Introduction 

O 



The work of Wyner [1J led to the development of the notion of secrecy capacity, which quan- 
tifies the maximum rate at which a transmitter can reliably send a secret message to a receiver, 



> 

•i-H 

X 

without an eavesdropper being able to decode it. More recently, researchers have considered 
secrecy for the two-user broadcast channel, where each receiver acts as an eavesdropper for the 
independent message transmitted to the other. This problem was addressed in [2J, where inner 
and outer bounds for the secrecy capacity region were established. Further work in Q studied 
the multiple-input single-output (MISO) Gaussian case, and [0 considered the general MEMO 
Gaussian case. It was shown in [4] that, under an input covariance constraint, both confidential 
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messages can be simultaneously communicated at their respective maximum secrecy rates, where 
the achievablity is obtained using secret dirty-paper coding (S-DPC). However, under an average 
power constraint, a computable secrecy capacity expression for the general MIMO case has not 
yet been derived. In principle, the secrecy capacity for this case could be found by an exhaustive 
search over the set of all input covariance matrices that satisfy the average power constraint 
flU. Clearly, the complexity associated with such a search and the implementation of dirty-paper 
encoding and decoding make such an approach prohibitive except for very simple scenarios, and 
motivates the study of simpler techniques based on linear precoding. 

While low-complexity linear transmission techniques have been extensively investigated for 
the broadcast channel (BC) without secrecy constraints, e.g., 0-Q, there has been relatively 
little work on considering secrecy in the design of linear precoders for the BC case. In JH, 
we considered linear precoders for the MIMO Gaussian broadcast channel with confidential 
messages based on the generalized singular value decomposition (GSVD) flU, ifTOl . It was shown 
numerically in [[8l that, with an optimal allocation of power for the GSVD-based precoder, the 
achievable secrecy rate is very close to the secrecy capacity region. 

In this paper, we show that for a two-user MIMO Gaussian BC with arbitrary numbers of 
antennas at each node and under an input covariance constraint, linear precoding is optimal and 
achieves the same secrecy rate region as S-DPC for certain input covariance constraints, and 
we derive an expression for the optimal precoders in these scenarios. We then use this result 
to develop a sub-optimal closed-form algorithm for calculating linear precoders for the case of 
average power constraints. Our numerical results indicate that the secrecy rate region achieved 
by this algorithm is close to that obtained by the optimal S-DPC approach with a search over 
all suitable input covariance matrices. 

In Section ||l[ we describe the model for the MIMO Gaussian broadcast channel with confi- 



dential messages and the optimal S-DPC scheme, proposed in [|4j|. In Section III we consider a 
general MIMO broadcast channel under a matrix covariance constraint, we derive the conditions 
under which linear precoding is optimal and achieves the same secrecy rate region as S-DPC, 
and we find the corresponding optimal precoders. We then present our sub-optimal algorithm for 



designing linear precoders for the case of an average power constraint in Section IV followed 
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by numerical examples in Section |V| Section [VI] concludes the paper. 



Notation: Vector- valued random variables are written with non-boldface uppercase letters (e.g., 
X), while the corresponding non-boldface lowercase letter (x) denotes a specific realization of 
the random variable. Scalar variables are written with non-boldface (lowercase or uppercase) 
letters. The Hermitian (i.e., conjugate) transpose is denoted by (.) H , the matrix trace by Tr(.), 
and I indicates an identity matrix. The inequality A y B (A y B) means that A — B is 
Hermitian positive (semi-)definite. Mutual information between the random variables A and B 
is denoted by I (A; B), E is the expectation operator, and CJ\f(0,a 2 ) represents the complex 
circularly symmetric Gaussian distribution with zero mean and variance a 2 . 

II. Broadcast Channel and S-DPC 

We consider a two-receiver multiple-antenna Gaussian broadcast channel with confidential 
messages, where the transmitter, receiver 1 and receiver 2 possess n t , m l5 and m 2 antennas, 
respectively. The transmitter has two independent confidential messages, W\ and Wi, where W\ 
is intended for receiver 1 but needs to be kept secret from receiver 2, and W 2 is intended for 
receiver 2 but needs to be kept secret from receiver 1 H|. 

The signals at each receiver can be written as: 

yi = Hx + zi 

(1) 

y 2 = Gx + z 2 

where x is the n t x 1 transmitted signal, and z, G C miXl is white Gaussian noise at receiver i 
with independent and identically distributed entries drawn from CJ\f(0, 1). The channel matrices 
H e C miXnt and G e Q"i 2 xn t are assumec i t be unrelated to each other, and known at all three 
nodes. The transmitted signal is subject to an average power constraint when 

Tt(E{XX h }) = Tr(Q) < P t (2) 

for some scalar P t , or it is subject to a matrix power constraint when [0, [fTTTl : 

E{XX H } = Q 1 S (3) 
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where Q is the transmit covariance matrix, and S >z 0. Compared with the average power 
constraint, ([3]) is rather precise and inflexible, although for example it does allow for the 
incorporation of per-antenna power constraints as a special case. 

It was shown in [2j that for any jointly distributed (Vx,V 2 ,X) such that (Vi,V 2 ) — > X — » 
(Yi, Y 2 ) forms a Markov chain and the power constraint over X is satisfied, the secrecy rate pair 
(Ri,R 2 ) given by 

R 1 = I(V 1 ;Y 1 )-I(V 1 ;V 2 ,Y 2 ) 

(4) 
R 2 = I(V 2 ;Y 2 )-I(V 2 ;V 1 ,Y 1 ) 

is achievable for the MIMO Gaussian broadcast channel given by ([TJ, where the auxiliary 
variables V\ and V 2 represent the precoding signals for the confidential messages W\ and W 2 , 
respectively 10. In O, the achievablity of the rate pair (|4]) was proved. 

Liu et al. [ffl analyzed the above secret communication problem under the matrix power- 
covariance constraint §3§. They showed that the secrecy capacity region C S (H, G,S) is rectan- 
gular. This interesting result implies that under the matrix power constraint, both confidential 
messages W\ and W 2 can be simultaneously transmitted at their respective maximal secrecy 
rates, as if over two separate MIMO Gaussian wiretap channels. To prove this result, Liu et al. 
showed that the secrecy capacity of the MIMO Gaussian wiretap channel can also be achieved 
via a coding scheme that uses artificial noise and random binning [4, Theorem 2]. 

Under the matrix power constraint ([3]), the achievablity of the optimal corner point (R[, R^) 
given by [4, Theorem 1] 

R*,= max log |HKiH H + l| -log|GK i G i/ + l| 



m = log 



GSG^+I 



(5) 



+ RI 



USU H + 1 

is obtained using dirty -paper coding based on double binning, or as referred to in |@]|, secret 
dirty paper coding (S-DPC). More precisely, let K^ y maximize d5|), and let 

V 1 = U 1 + ¥U 2 V 2 = U 2 X = U 1 + U 2 , (6) 
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where U\ and U 2 are two independent Gaussian vectors with zero means and covariance matrices 
K* t and S— K*, respectively, and the precoding matrix F is defined as F = K^H^ (HK;H fl + i) H. 
One can easily confirm the achievablity of the corner point {R[,R* 2 ) by evaluating ^ for the 
above random variables and noting that in ([TJ), X = U\ + Ui. Note that under the matrix power 
constraint S, the input covariance matrix that achieves the corner point in the secrecy capacity 
region satisfies Q = S (J. 

The matrix K t that maximizes ([5]) is given by flU, [fTTTl 



K? 



S^C 



(Cfd)- 1 



C"S^ 



(7) 



Q Q 

where C = [Ci C2] is an invertiblej generalized eigenvector matrix of the pencil 



(S^H^HS^ + I , S^G^GS^ + i) 



(8) 



satisfying 02 






S5H"HS2 +1 

S^G H GS^+I 



C = A 
C = I, 



(9) 



where A = diag{Ai, ..., X nt } >- contains the generalized eigenvalues sorted without loss of 
generality such that 

Ai > ... > \ b > 1 > \ b+ i > ... > A ni > . 

The quantity b denotes the number of generalized eigenvalues greater than one (0 < b < n t ), 
and defines the following matrix partitions: 



Ai 

A 2 



[Ci c 



2J , 



(10) 



where A x = diag{Ai, ..., At}, A 2 = diag{A(, + i, ..., A nt }, Ci contains the b generalized eigenvec- 
tors corresponding to Ai and C 2 the (n t — b) generalized eigenvectors corresponding to A 2 . 



Note that C is invertible since both components of the pencil \S\ are positive definite. 
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Now, by applying ^ in ^, the corner rate pair (R* 1: R* 2 ) can be calculated as ([4, Theorem 3]) 

^?i = log | Ai| 

(11) 
i^ = -log|A 2 | . 

For the average power constraint in Q, there is no computable secrecy capacity expression 
for the general MIMO case. In principle the secrecy capacity region for the average power 
constraint, C S (H, G,P t ), could be found through an exhaustive search over all suitable matrix 
power constraints [|4l. [[T3l Lemma 1]: 

C s (H,G,P t ) = |J C S (H,G,S). (12) 

SbO,Tr(S)<P t 



For any given semidefinite S, C S (H, G, S) can be computed as given by ( 1 1 ). Then, the secrecy 
capacity region C S (H, G, P t ) is the convex hull of all of the obtained corner points using (11). 
The complexity associated with such a search, as well as that required to implement dirty- 
paper encoding and decoding, are the main drawbacks of using S-DPC to find the secrecy 
capacity region C S (H,G, P t ) for the average power constraint. This makes linear precoding 
(beamforming) techniques an attractive alternative because of their simplicity. To address the 
performance achievable with linear precoding, we first describe the conditions under which linear 
precoding is optimal in attaining the same secrecy rate region that is achievable via S-DPC, when 
the broadcast channel is under an input covariance constraint. In particular, in the next section 
we show that this equivalence holds for matrix power constraints that satisfy a certain property, 



and we derive the linear precoders that achieve optimal performance. Section IV then uses these 



results to derive a sub-optimal algorithm for the case of the average power constraint. 

III. Optimality of Linear Precoding for BC Secrecy 
In this section we answer the following questions: 

(a) For a given general MEMO Gaussian BC described by ([T]), where each node has an arbitrary 
number of antennas and the channel input is under the covariance constraint ([3]), is there 
any S y for which linear precoding can attain the secrecy capacity region? 

(b) If yes, how can such S be described? 
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(c) For such S, what is the optimal linear precoder that allows the rectangular S-DPC capacity 



region given by ( 1 1 ) to be achieved? 
(d) If S does not satisfy the condition for optimal linear precoding in (a), what is the worst-case 

loss in secrecy capacity incurred by using the linear precoding approach described in (b) 

anyway? 
To begin, we give the following theorem as an answer to questions (a) and (b) above. 

Theorem 1. Suppose the matrix power constraint S ^ on the input covariance Q in ([3]) leads 
to generalized eigenvectors in (J9j) that satisfy span{Ci} _L span{C 2 }, i.e. CfC 2 = 0. Then the 
secrecy capacity region C S (H, G, S) can be achieved with X = V\ + V 2 , where V\ and V 2 are 
independent Gaussian precoders respectively corresponding to W\ and W 2 , with zero means and 
covariance matrices K^ and S — K.* t , with K^ defined in (|7]). 

Proof: Recall that for any S >z 0, the secrecy capacity region C S (H, G, S) is rectangular, 
so we only need to show that when C^Ci = 0, the linear precoders V\ and V 2 characterized 
in this theorem are capable of achieving the corner point (i?*, i? 2 ) given by (11 ). From Q, the 
achievable secrecy rate R\ is given by 

R x = /(Vi; Fi) - I(Vx; V 2 ) - J(Vi; Y 2 \ V 2 ) = I(Vy, Y x ) - 7(y l5 Y 2 \V 2 ) (13) 
= I(V X ; U(Vx + V 2 ) + Z x ) - I(V X ; G(V X + V 2 ) + Z 2 |V 2 ) 

= J(y i; H(V X + V 2 ) + Z X ) - I{V X - GV X + Z 2 ) (14) 

= log|A 1 | = i2* , (15) 



where (13) and the second part of (14) come from the fact that V x and V 2 are independent. 
Equation (15) is proved in Appendix A. One can similarly show that R 2 = i? 2 = — log | A 2 | is 
achievable to complete the proof. ■ 

Theorem [T] shows that the secrecy capacity region corresponding to any S with orthogonal 
Ci and C 2 can be achieved using either linear independent precoders V\ and V 2 , as defined in 
Theorem [T] or using the S-DPC approach, as given by ([6]). The next theorem expands on the 
answer to question (b) above, and also addresses (c). First however we present the following 
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lemma which holds for any S y 0. 



Lemma 1. For a given BC under the matrix power constraint ([3]), for any S >z we have 
rank(Ci) < m, where m is the number of positive eigenvalues of the matrix H H H — G^G. 

Proof: Please see Appendix B. ■ 

The following theorem presents a more specific condition on S that results in generalized 
eigenvectors that satisfy span{Ci} _L span{C 2 }. 

Theorem 2. For any S >z 0, the generalized eigenvectors Ci and C2 in (|9]) are orthogonal iff 
there exists a matrix T G C ntXnt such that S = TT^ and T simultaneously block diagonalizes 
H^H and G H G: 



rH^HT 



K H i Q 

Q K H2 



rG w GT 



K G1 
Q K G2 



(16) 



where the m x m matrices K H i h and K G i h satisfy K H i h K G1 and K H2 d K G2 . 
Proof: The proof begins by noting that if S = TT H , then the pencil in d8J) and 



(T^H^HT + I , T F G H GT + I) 

have exactly the same generalized eigenvalue matrix A, and thus the same secrecy capacity 
regions. The remainder of the proof can be found in Appendix C. ■ 

While algorithms exist to find T that jointly block diagonalizes H H H and G^G (see for 
example llT4l and references therein), as mentioned in Appendix C only those T that lead to 
Khi b K G1 and K H2 d> K G2 are acceptable. Later, we will demonstrate that for any BC there 
are an infinite number of matrix constraints S that can achieve such a block diagonalization and 
hence allow for an optimal linear precoding solution. 

To conclude this section, we now answer question (d) posed above. Define the projection 
matrices P G . = Cj(Cf f Cj) _1 C,f / and P Cj = I — P Gj , and note that in general, equation a7| 
is equivalent to K* = S2P Cl S5 and S - K* = S5P Ci S2. When Cf C 2 = 0, the optimal 
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covariance matrices for V\ and V 2 also satisfy 



K* = S^Pi S^ (17) 



S-K* = S5P C2 S2 . (18) 

The following theorem explains the loss in secrecy that results when linear precoding with these 
covariances is used for a matrix constraint S that does not satisfy C^C 2 = 0. 

Theorem 3. Assume a linear precoding scheme X = V\ + V 2 for independent Gaussian precoders 
V\ and V2 with zero means and covariance matrices K 4 = S^P^S 2 and S — K 4 = S2P C2 S5, 
respectively. Also define N = (C^Pq i C 2 ) C^P^P^^i- The loss in secrecy capacity that 
results from using this approach in the two-user BC is at most log |l + N^N| for each user. In 
particular, the following secrecy rate pair is achievable: 

R x = max(0, R{ - log ll + N^Nl) = max(0, log |Ai| - log ll + N H N\) 

(19) 
R 2 = max(0, R* 2 - log |l + N H N\) = max(0, -log |A 2 | - log |l + N H N\) . 

Proof: See Appendix D. ■ 

Remark 1. Note that if Ci and C 2 are orthogonal, then N = and Ri = R*, i = 1,2, is 
achievable, as discussed in Theorem 1. 

IV. Sub-Optimal Solutions Under an Average Power Constraint 

So far we have shown that if the broadcast channel ([T]) is under the matrix power constraint 
S Q, then linear precoding as defined by Theorem [T] is an optimal solution when S satisfies the 
condition described in Theorem [2} In the following we propose a suboptimal closed-form linear 
precoding scheme for the general MIMO Gaussian BC under the average power constraint ([2]), 
where as mentioned earlier there exists no optimal closed-form solution that characterizes the 
secrecy capacity region. We begin with some preliminary results, then we develop the algorithm 
for the general MEMO case, and finally we present an alternative algorithm specifically for the 
MISO case since it offers additional insight. 



April 18, 2013 DRAFT 



10 



A. Preliminary Results 

Remark 2. Suppose that the input covariance matrix Q leads to a point on the Pareto boundary 
of the secrecy capacity region given by (12) under the average power constraint ([2]). Then 
Tr(Q) = P t and Q cannot have any component in the nullspace of H H H + G^G, and thus 
C.(H, G, P t ) = C s (U eq , G eq , P t ), where H eq 
vectors corresponding to the non-zero singular values of H H H + G^G. 



HU P , G eq = GU P and U p contains the singular 



According to Remark [2} we can assume without loss of generality that H H H + G H G is full- 
rank; otherwise, we could replace H, G with H eq , G eq and have an equivalent problem where 
H^H e q + G^G ei j is full-rank and the secrecy capacity region is the same (in such a case, n t 
would then represent the number of transmitted data streams rather than the number of antennas). 
With this result, we have the following lemma. 



Lemma 2. Define 



W 



tH 



H^H + G^G 



Hi 



(20) 



Then WH H HW and WG fl GW commute and hence share the same set of eigenvectors: 

WH^HW = $ W S X $^ 
WG ff GW = $ W S 2 $^ , 



(21) 



where <fr w is the (unitary) matrix of eigenvectors and S x y 0, S 2 ^0 the corresponding 
eigenvalues. 

Proof: See Appendix E. ■ 

Without loss of generality, we assume that the columns of <fr w are sorted such that the first 

p diagonal elements of Si are greater than the first p diagonal elements of S 2 , and the last 

n t — p diagonal elements of Si are less than or equal to those of S 2 . Recall from Lemma [T] 

that < p < m, where m is the number of positive eigenvalues of H H H — G H G. Thus, 



Si 











'1/5 



^2p 









J 2p 



(22) 
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where S ip is p x p, T, ip is (ra t - p) x (n t - p), £ 1p >- £ 2p and S 1? ^ £ 2p -. 
Now define 






w J-w 



W$ w P $"w 



W$ W P2 



(23) 
(24) 



where W and <fr w are given in Lemma [2] and P ^ is any block-diagonal matrix partitioned 



in the same way as Si and £ 2 . With these definitions, we see from (21 ) that T^H HT W and 



T^G^GT W are block diagonal. Thus, from Theorem [2| a BC with the matrix power constraint 
S w = T^T W leads to a matrix pencil (sIh^HsI + 1 , sIg^GsI + 1) with generalized 
eigenvectors C w = [C iw C 2w ] that satisfy Cf^C 2w = 0, where Ci w , C 2w correspond to 
generalized eigenvalues that are larger or less-than-or-equal-to one, respectively. 

Remark 3. Since the above result holds for any block-diagonal P y with appropriate 
dimensions, then for every BC there are an infinite number of matrix power constraints S w 
that achieve a block diagonalization and hence allow for an optimal linear precoding solution. 

In the following, we restrict our attention to diagonal rather than block-diagonal matrices P, 
for which a closed form solution can be derived. From Theorem [T| we have the following result. 

Lemma 3. For any diagonal P y 0, the secrecy capacity of the broadcast channel in ([T]) under the 



matrix power constraint S w = T W T^ defined in (20)-(24) can be obtained by linear precoding. 
In particular, 



X = W$, 



vi 



v 1 + v 2 



(25) 



where V[ G C p and V^ E C nt p are independent Gaussian random vectors with zero means and 
covariance matrices Pi and P 2 such that 



Pi 
P 2 



(26) 



and as before Vi, V 2 represent independently encoded Gaussian codebook symbols corresponding 
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to the confidential messages W\ and W 2 , with zero means and covariances K^ w and S w — K^ w 
respectively given by 



K 



tw 



K 



tw 



W$ Y 



W$ v 



Pi 




















p 2 



$"w 



*"w. 



(27) 



(28) 



Proof: The matrix T w simultaneously block diagonalizes H H H and G H G, so by Theo- 
rems [T] and [2] we know that linear precoding can achieve the secrecy capacity region. The proof 



is completed in Appendix F by showing the equality in (25), and showing that (27) corresponds 
to the optimal covariance in Q. ■ 



From the proof in Appendix F and (21 )-(24), we see that under the matrix power constraint 



S w given by (23 ) with diagonal P, the general BC is transformed to an equivalent BC with a set 
of parallel independent subchannels between the transmitter and the receivers, and it suffices for 
the transmitter to use independent Gaussian codebooks across these subchannels. In particular, 
the diagonal entries of Pi and P 2 represent the power assigned to these independent subchannels 



prior to application of the precoder W<fr w in (|24|p| From (25), the signals at the two receivers 
are given by 



y-i 



HW$, 



TiEi 



Zl 



Ty 



Y'2 



Sipv; 

s 2 pv; 
s 2p -v^ 



Zl 



+ Zi 



+ z 2 



2 Note that the matrices Pi, P2 do not represent the actual transmitted power, since the columns of W# w are not unit-norm. 
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where r 1; T 2 are unitary. The confidential message for receiver 1 is thus transmitted with power 
loading Pi over those subchannels which are degraded for receiver 2 (S ip y S 2p ), while 
receiver 2's confidential message has power loading P 2 over subchannels which are degraded 
for receiver 1 (S 2p - >- Sip). Any subchannels for which the diagonal elements of S 2p - are equal 
to those of Sip are useless from the viewpoint of secret communication, but could be used to 
send common non-confidential messages. 

From Theorem [Tj the rectangular secrecy capacity region of the MIMO Gaussian BC ([T]) under 
the matrix power constraint S w (|23~|) is defined by the corner points 



i^(Pi) = log |A lw | = log |I + SipPil - log |I + S 2 pP!| 

(29) 

R* 2 (P 2 ) = - log |A 2w | = log |I + S 2p -P 2 | - log |I + SipP 2 | , 



where A iw is given by ( 8 1 ) in Appendix F. Note that we have explicitly written R\ as a function 
of the diagonal matrix P 1 y to emphasize that Pi contains the only parameters that can be 
optimized for R^. More precisely, since for a given matrix power constraint S w , S ip and S 2p 



are channel dependent and thus fixed, as shown in (21 )-(22). A similar description is also true 



for R*. 



B. Algorithm for the MIMO Case Under the Average Power Constraint 

Here we propose our sub-optimal closed form solution based on linear precoding for the 
broadcast channel under the average power constraint ([2]). The goal is to find the diagonal 



matrix P in (23) that maximizes R* in (29) for a given allocation of the transmit power to 



message W{, and that satisfies the average power constrainjj 

Tr{E{XX H }) = Tr(S w ) = Tr(W$ w P$^W) 

= Tr($^W 2 $ w P) = Tr ($^(H H H + G // G)- 1 $ w P) = P t . (30) 



3 Note that since we want to characterize the achievable secrecy rate points on the Pareto boundary, we use an equality 
constraint on the total power P t in j30b. 
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Noting that <fr w can be written as <fr w = [<&i w 3?2w], where <&i w is a n t X p submatrix 



corresponding to the eigenvalues in Si p , (30) can be rewritten as 

Tr(£{XX^}) = Tr ($£(H H H + G ff G) ^ w P) 

= Tr ($f w (H"H + G H G) ^ Pj + Tr ($f w (H H H + G H G)- 1 $ 2w P 2 ) 
= Tr(A 1 P 1 ) + Tr(A 2 P 2 ) = P t (31) 

where we defined positive definite matrices Aj = $^(H H H + G^G) _1 $j W , i — 1,2. 

Our sub-optimal closed-form solution for the BC under the average power constraint ([2]) is 
not optimal, since instead of doing an exhaustive search over all S y with Tr(S) = P t as 
indicated in ( p"2] ), we will only consider specific S matrices of the form given for S w in ( |2~3| ) 
with diagonal P. Since i?*(Pj) is only a function of P i5 i£*(Pi) and i££(P 2 ) can be optimized 
separately for any power fraction a (0 < a < 1) under the constraints Tr (A]Pi) = a? t and 
Tr (A 2 P 2 ) = (1 - a)P t , respectively. 

Theorem 4. For any a, < a < 1, the diagonal elements of the optimal P^ and P2 are given 

by 

* ( n — {(Tipi + o- 2p i) + VA^V ~~ °"2 P i) 2 + 4(<Ti p i - cr 2p j)cripiCr 2p j/(/iiaij) \ 
\ z cr lp jO" 2p j / 

* _ /„ — (o"ipi + a 2 pi) + y/((T2pi - v\pi) 2 + ^(cr 2 pi - Vipi)v 2 pi(Jipi/ (p 2 a 2i ) \ 

Vli — maX U ' o i (33) 

\ 2 a lpi a 2 pi J 

where a^, a 2pi ,, and an are the 2 th diagonal elements of H lp , £ 2p , and Ai, respectively, where 
< i < p. Also a\pi, a 2 pi, and a 2 i are the i th diagonal elements of Sip, £ 2/5 , and A 2 , respectively, 
where < i < [n t — p). The Lagrange parameters \X\ > and /i 2 > are chosen to satisfy the 
average power constraints Tr (A1P1) = aP t and Tr (A 2 P 2 ) = (1 — a)P t , respectively. 

Proof: We want to optimize diagonal matrices Pi and P 2 so that the secrecy rates R\(P\) 



and Rl(P 2 ), given by (29), are maximized for a given a, < a < 1. Since i?*(Pj) only 



depends on Pj, the two terms in (29) can be maximized independently. We show the result 



for i — 1; the procedure for i = 2 is identical. From (29), the Lagrangian associated with 
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max Tr(Al p l)=aPt i?i(Pi) is 



£ = log |I + S lp Pi| - |I + H 2p P 1 \ - //iTrCAiPi) 

= ^[log(l + aipipii) - log(l + a 2p iPii)] - H\ 22 a uPii > 



(34) 



where ^x > is the Lagrange multiplier. Since S lp >- S 2p , Eq. (34) represents a convex 



optimization problem. The optimal P^ with diagonal elements given by (32) is simply obtained 



by applying the KKT conditions to (34). 



Corollary 1. For any a, < a < 1, let Rl(a) and R 2 (a) represent the corner points given by 



< |29| ) for the optimal P^ and P^, given by ( |32| ) and < |33[ ). The achievable secrecy rate region of 
the above approach under the average power constraint ([2]) is the convex hull of all obtained 
corner points and is given by 



TZ s (H,G,P t )= |J (R{(a) , R* 2 (a)) . 



(35) 



0<a<l 



It is interesting to note that, unlike the conventional broadcast channel without secrecy con- 
straints where uniform power allocation is optimal in maximizing the sum-rate in the high SNR 
regime lfT6ll . the high SNR power allocation for the BC with confidential messages is a special 
form of waterfilling as described in the following lemma. 



Lemma 4. For high SNR (P t — ¥ oo), the asymptotic optimal power allocations given by (32) 



(33) are 



Pu 



P2i 



1 



1 



/ilOij \(J2pi 0\ 



I" 



(36) 



1 



1 



1 



/I2CI2; V "!^ °"2pi 



(37) 



Pi 



Proof: To show (36) we note that /ii — )• when P t — > 00. Thus ( [32] ) can be written as 
■\/Hp\pi - cr 2pi )<7 lpi a 2p i/(lJ,ia u } _ 



2 <3\pi<J 2 pi 



1 cripi — a 2p i 



1 



1 



1 



\l\Q,\i \0 2p i 0\ p i 



(37) is proved similarly. 
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It is also worth noting that the solution in p2|)-(|3~3~|) approaches the standard point-to-point 



MIMO waterfilling solution when one of the channels is dominant. For example, let G — > 0. 
We will show that the optimal input covariance simplifies to the waterfilling solution for H, 
given by Q* P = <fr# (-1 — A // 1 ) + <fr^, where unitary <&# and diagonal A# are obtained from 
the eigenvalue decomposition H H H = <& H A H <&x- The capacity of the point-to-point MIMO 
Gaussian link is 



C = log I + HC&H 



H\ 



log 



I 



A,/ ( -I 



A^ 1 



(38) 



When G ->• 0, we note from (20) and (21) that $ w ->■ <& H , S x = E 1/0 ->• I, S 2 ->■ 0, and 



Pi = P. Consequently, i?2 — > and iZ* -» log|I + P*|, where P* is a diagonal matrix with 



diagonal elements given by (32). The average power constraint in (31 ) becomes Tr (AP) = P t , 
where A = $^(H H H + G^G) x $ w -^ A^ 1 when G -^ 0. Thus the 2 th diagonal element of 



A converges to the i th diagonal element of A^ 1 . Starting from (32) and applying L'Hopital's 
rule, when G — > and hence a 2p i — V 0, we have P* —¥ A#(-I — A^) + , and consequently, 



lim R\ = log |I + P* | = log 

G->0 



I + Anf-I-A^ 1 



C. Alternative Approach for the MISO Case 

Here we focus on the BC in ([TJ) for the MISO case under an average power constraint, where 
both receivers have a single antenna and the transmitter has n t > 2 antennas: 

2/i = h^x + z\ 



m 



x + z 2 , 



where the channels are represented by the nt x 1 vectors h and g. The MISO case is the only BC 
scenario whose secrecy capacity region under the average power constraint ([2]) is characterized 
in closed-form. In particular, it was shown in Q that 



C s (h,g,P t )= |J {C x {a),C 2 {a)) 



(39) 



0<a<l 
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where (Ci(a) , C 2 (a)) is the secrecy rate pair on the Pareto boundary of the secrecy capacity 
region for the power fraction a, < a < 1, where power aP t is allocated to receiver l's message 
and (1 — a)P t is allocated to receiver 2's message. Furthermore, we have 

Ci(a) = log7i(a) 

(40) 

C 2 {a) = log 72 (a) , 

where 

1 + aPt ef hh^ei 
lM) ~ 1 + aP t ef gg^ ei ' 

ei is the unit length principal generalized eigenvector of (I + P t hh H ; I + P t gg H ), 72(a) is the 
largest generalized eigenvalue of 

1 l + «P t |efg| SS ' l + aP t \efh\ ' 

and e 2 denotes the unit length generalized eigenvector corresponding to 72(a)- Note that the 



achievablity of (40) is still based on S-DPC. 



While we could have just used the results of Section IV-B for the MISO case, we will see 
that the advantage of considering a different approach here is that we obtain a more succinct 
expression for the achievable secrecy rate region for linear precoding, and we are able to quantify 
the loss in secrecy rate incurred by linear precoding under the average power constraint compared 
with (Ci(a),C 2 (a)). This was not possible in the MIMO case. 



Referring to ([6]), it was shown in [[3l that for the secrecy rate pair given by (40), U\ and 



U 2 have covariance matrices aP t e,ie^ and (1 — a)P t e 2 e 2 Y , respectively. Thus, the specific input 



covariance matrix that attains (40) is given by 



S Q = aP t e^ + (1 - a)P t e 2 e^ , (41) 

where Tt(Sq) = P t and rank(Sg) = 2. Equivalently, one can say that under the matrix power 
constraint Sq, the corner point of the corresponding rectangular secrecy capacity region is given 



by (40). The union of these corner points constructs the Pareto boundary of the secrecy capacity 



region under the average power constraint, where any point on the boundary is given by d40 
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for a different a and is achieved under the matrix power constraint Sq given by (41 ). 

Using the above fact, we now present a different linear precoding scheme as an alternative to 
Corollary [T] for the MISO BC under the average power constraint Q. 

Corollary 2. Using the linear precoding scheme proposed in Theorem [3] for the MISO BC under 
an average power constraint, the following secrecy rate region is achievable: 

K.(h,g,P t ) = |J (R l (a),R 2 (a)) , 



0<a<l 



where 



R^a) = max^a) - log (l + (cf P^c 2 )- 2 |cf P£,P£c 2 | 2 ) , 0) 
R 2 (a) = max(C 2 (a) - log (l + (cf P^c,)" 2 |ef P^P^j 2 ) , 0) , 



Ci(a) and C 2 (a) are given by (40), 



(42) 



c-i 



c 2 



eftS- 1 



Q 



r H 



S^ ei 



ei 



So'fi 



(43) 
(44) 



^(S^+gg^f, 

and where fi is the unit length principal generalized eigenvector of (I + Ptgg H ; I + Pthh H ). 

Proof: From Remark [2} and by noting that for any MISO BC, hh^ + gg H has at most 2 
non-zero eigenvalues, any MISO BC can be modeled with a scenario involving just two transmit 
antennas. Thus, without loss of generality, we assume that n t = 2. From Theorem |3j we only 
need to characterize Ci and c 2 , where Ci (c 2 ) is the generalized eigenvector of the pencil 



s|hh"s| + I , S|: 



M 



SA + I 



(45) 



corresponding to the generalized eigenvalue larger (less) than 1, Ai (A 2 ). 
From (|6]) and Q, the covariance matrix of U\ can be rewritten as 



K? 



i 
SA 



[ci c 2 ] 



(Cfc!)" 1 








cic^rsi 



C2„ n H G2 



(46) 
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Fig. 1. Secrecy capacity region of S-DPC together with secrecy rate region for linear precoding with P t = 12, tit = mi 
mi = 2. 



Comparing (46) with the covariance matrix of U\ reported in 0, we have aP^ie^ = r^ SqCi c^Sq. 



This results irQ 



ci 

l c i| 



aP t S 2 e 1 



(47) 



On the other hand, from the definition of Ci and c 2 (see (|9])-([T0|) for example), we have 



Ci c 2 



Ci c 2 



H 



H 



S Jhh"S J + 1 



S^gg^S| + 1 



Ci c 2 



ci c 2 



Ai 
A 2 



7i(a) 
^\a) 



(48) 



where 71(a) and 72(a) are defined after (40), and the fact that Ai = 71(a) and A 2 



72(a) 



comes from the argument after (41 ) and by comparing (11 ) and (40). Substituting (47) in (48), 
after some simple calculations, c x can be explicitly written as in ( |43j ). Recalling that Ci is the 



principal generalized eigenvector of (45) and c 2 , which corresponds to the smallest generalized 



4 Note that multiplication by a factor exp(j9) is required for a precise equaltiy, but since this term disappears in the final 
result, we simply ignore it. 
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eigenvalue of the pencil ([45]), is the principal generalized eigenvector of the pencil 



(s|g g "sj + 1 , sjhh^sj + i) 



we obtain (44). The proof is completed by using (43) and (44) in (19) 



V. Numerical Results 




t* 1 1.5- 



1.5 2 2.5 3 

fi r (bps/Hz) 



3.5 



Fig. 2. Secrecy capacity region of S-DPC together with secrecy rate region for linear precodings in Cor. 1 and Cor.2, with 
Pt = 10, n t = 2,m 1 =m 2 = 1. 



In this section, we provide numerical examples to illustrate the achievable secrecy rate region 
of the MEMO Gaussian BC under the average power constraint ([2]). In the first example, we have 
P t = 12, H = [0.3 2.5; 2.2 1.8] and G = [1.3 1.2; 1.5 3.9], which is identical to the case studied 
in [|4l Fig. 3 (d)]. Fig. [T] compares the achievable secrecy rate region of the proposed linear 
precoding scheme in Section IV-A with the secrecy capacity region obtained by the optimal S- 
DPC approach together with an exhaustive search over suitable matrix constraints, as described 
in Section II. We see that in this example, the performance of the proposed linear precoding 
approach is essentially identical to that of the optimal S-DPC scheme. 

In the next example, we study the MISO BC for P t = 10. Fig. [2] shows the average secrecy 
rate regions for S-DPC and the suboptimal linear precoding algorithms described in Corollary 1 



April 18, 2013 



DRAFT 



21 

and 2. This plot is based on an average of over 30000 channel realizations, where the channel 
coefficients were generated as independent £/V(0, 1) random variables. We see that Corollary 2 
provides near optimal performance when a — ¥ or a — > 1, while Corollary 1 is better for 
in-between values of a. The degradation of using linear precoding with Corollary 1 is never 
above 15% for any a. 

VI. Conclusions 

We have shown that for a two-user Gaussian BC with an arbitrary number of antennas at 
each node, when the channel input is under the matrix power constraint, linear precoding is 
optimal and achieves the secrecy capacity region attained by the optimal S-DPC approach if the 
matrix constraint satisfies a specific condition. We characterized the form of the linear precoding 
that achieves the secrecy capacity region in such cases, and we quantified the maximum loss 
in secrecy rate that occurs if the matrix power constraint does not satisfy the given condition. 
Based on these observations, we then formulated a sub-optimal approach for the general MEMO 
scenario based on linear precoding for the case of an average power constraint, for which no 
known characterization of the secrecy capacity region exists. We also studied the MISO case in 
detail. Numerical results indicate that the proposed linear precoding approaches yield secrecy 
rate regions that are close to the secrecy capacity achieved by S-DPC. 



Appendix A 



Proof ofEq. (15) 



From (14), we have 



R x = I(V i; H(K + V 2 ) + Z x ) - I(V i; GVt + Z 2 ) 

= log |I + SH^Hl - log |I + (S - K*)H H H| - log ll + K*G H G| 



(49) 



The covariance K* t , given by (|7]), can be rewritten as 



k: 



S*[Ci C 2 ] 



(Cfd)" 1 







s^ 
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S^ d(Cf CO^Cf S3 = S^ P Cl s^ , 



(50) 



where P Cl = Ci(CfCi) 1 Cf / is the projection matrix onto the column space of Ci. Moreover, 
let Pq = I — Pd be the projection onto the space orthogonal to Ci. Consequently, we have 



S-K* = S-S2P Cl S 



S2P^ S2 =S2Pn,S2 



s^c 



Ci >J " ~ >J " x c 2 

Q Q 

Q (Cfd)- 1 



C"S^ 



(51) 

(52) 



where in (51 ), P^ = Pc 2 comes from the fact that span{Ci} JL span{C 2 }, and C = [Ci C 2 ] 
is full-rank. 

Following the same steps as in the proof of |[P3l Lemma 2] or 0] App. B], we can convert 
the case when S ^ 0, |S| = 0, to the case where S >- with the same secrecy capacity region. 
From (|9]) and ( flO] ) we have 



u H u = s 1 / 2 



c 



-H 



Ai 
A 2 



c 



-1 



I 



g-1/2 



(53) 



G H G = S" 1 / 2 [C-^C- 1 - I] S- 1 / 2 



Using (52) and (53), we have: 



|I + fS - K?)H H H 



I + S^C 



1 + 



C"- 



Q Q 

Q (CfC.)- 1 

Q Q 
Q (CfCa)- 1 

I Q 

(C 2 ¥ C 2 )- 1 A 2 

(CfC 2 )-^ 2 | = |(CfC 2 )-''-|A 



C 



11 



Ai 
A 2 



C 



-1 



Ai 
A 2 



-C^C 



g-1/2 

(54) 

(55) 
(56) 



where (54) comes from the fact that |I + AB| = |I + BA|. Finally, (55) holds since Cf C 2 = 
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and C^C is block diagonal. 
Similarly, one can show that 



|l + K*G H G| 



Kcfd 



and 



I+SH W H 



H-, 



c 



-H 



A t 
A 2 



i-i 



flri\-l| 



KcfCj 



KCfCa)- 1 ! 



fC"C 



|Ai| • |A 2 



|Ai| • |A 2 



(57) 



(58) 



Substituting (56), (57) and (58) in (49), we have R x = log |Ai| = R\, and this completes the 
proof. 



Appendix B 
Proof of LemmaQ] 

From (|9|)-([T0|), we know that rank(Ci) = b, where b represents number of generalized 
eigenvalues of the pencil ([8]) that are greater than 1. From @-([10]), we have 



Cf 

cf 



S5H"HS2 +1 



S^G^GSa +1 



C^Ax 
Ci = I. 



(59) 
(60) 



Subtracting (59) from (60), a straightforward computation yields 



Cf S3 [U H U - G H G] S^d = Ai - I y . 



(61) 



From (6lJ, we have Cf Ss [H^H - G H G] S^Ci >- 0, from which it follows that rank(Ci) = 
b < m. Similarly one can show that rank(C 2 ) < m', where C 2 corresponds to the generalized 
eigenvalues of the pencil ([8]) which are less than 1, and m! represents number of negative 
eigenvalues of H^H - G H G. 
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Appendix C 
Proof of Theorem [2] 

We want to characterize the matrices S >z Q for which 



S2H W HS2 + I , S^G^GS^ + I J (62) 

has generalized eigenvectors with orthogonal Ci and C 2 . For any positive semidefinite matrix 
S G C ntXnt , there exists a matrix T G C ntXnt such that S = TT^ 02l. More precisely, 
T = Sa^, where \l/ can be any n t x n 4 unitary matrix; thus T is not unique. 

Remark 4. Let the invertible matrix C and the diagonal matrix A respectively represent the 
generalized eigenvectors and eigenvalues of 

(T^H^HT + I , T^G^GT + I) , (63) 

so that 

C H [T^H^HT + 1] C = A 

(64) 

C [T^G^GT + I] C = I , 



where T = S^^> for a given unitary matrix \I/. By comparing (M) and (64), one can confirm 



that \I/C = C and A = A, where C and A are respectively the generalized eigenvectors and 
eigenvalues of (62), as given by ([9]). 

17 

Also note that, for any unitary \I/, C C = C^C. Thus, finding a S >z Q such that (62) has 



orthogonal Ci and C 2 (block diagonal C H C) is equivalent to finding a T, S = TT H , such 



TJ 

that (63) has orthogonal Ci and C 2 (block diagonal C C). 



The if part of Theorem [2] is easy to show. We want to show that if S = TT H and T 



simultaneously block diagonalizes H H H and G H G, as given by ( 16) such that K H i h. K G1 and 



Kh2 _^ Kq2, then C^C 2 = 0. From the definition of the generalized eigenvalue decomposition, 
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we have 



C"[T H H"HT + I]C = C" 



C"[T"G H GT + I]C = C H 



I + K H i Q 

I + K H2 

I + K G1 

I + K G2 



C 



Dx 
D 2 



c = 


I Q 




Q I 



(65) 



from which we have 



C = [Ci c 2 ] 



Cu 

c 22 



where the invertible matrix Cu and diagonal matrix Di are respectively the generalized eigen- 
vectors and eigenvalues of (I + K m ; I + K G1 ). Since K m h K G1 , then D x y I> which 
shows that Ci corresponds to generalized eigenvalues that are bigger than or equal to one. We 
have a similar definition for C 22 and diagonal matrix D 2 , corresponding to (I + Kh2 ; I + K G2 ). 
Finally, since C C is block diagonal, then C^C, where C is the generalized eigenvector matrix 
of ( [62] ), is block diagonal as well. This completes the if part of the theorem. 

In the following, we prove the only if part of Theorem [2| i.e., we show that if S >z 



results in (62) having orthogonal Ci and C 2 , then there must exist a square matrix T such that 
S = TT H and T^H^HT and T E G ff GT are simultaneously block diagonalized as in <[l6) 
with K H i h K G i and K H 2 _^ K G2 . 

Let SzG^GSz have the eigenvalue decomposition $ B S B $^, where $b is unitary and H B is 
apositive semidefinite diagonal matrix. Also let (1+T, b )~^b (l + SsH^HS^ J $ B (I+S B ) _ 5 
have the eigenvalue decomposition f^S^f, where $a is unitary and S^ is a positive definite 
diagonal matrix. One can easily confirm that [fT2l C = <&s(I + £ B )~2$ A and A = S^, where 



C and A are respectively the generalized eigenvectors and eigenvalues of (62). Also let C be 
ordered such that C = [Ci C 2 ], where Ci corresponds to the generalized eigenvalues bigger 
than (or equal to) 1. We have 



c H c = ^(i + s B )- 1 $4 



(66) 



From (66), C^C is block diagonal iff the unitary matrix <& A is block diagonal. Recall- 



April 18, 2013 



DRAFT 



20 



ing that $ A is the eigenvector matrix of (I + E^-^g f I + SzH^HSa j 0> B (I + E B )~5, a 
block diagonal $ A leads to (I + S B )~2$f (l + SsH^HSs J $ B (I + S B )"i and consequently 
$f S5H H HS5$ B must be block diagonal. Thus, if C H C is block diagonal, i.e., Cf C 2 = 0, 
there must exist a unitary matrix <fr B such that tfrfSaH^HSs^ and <fr B S2G H GS2<I? B are 
simultaneously block diagonal|j Letting T = S2<3? B results in (65), for which we must have 
Khi h Kgi and K H 2 if! K G 2, otherwise it contradicts the ordering of C = [Ci C 2 ]. This 
completes the proof. 



Appendix D 
Proof of Theorem [3] 



We need to prove that the secrecy rate pair (_R l5 i? 2 ) given by (19) is achievable. 



Remark 5. By applying the Schur Complement Lemma lfT2l on 



H 



C^C = [Ci C 2 ] [Ci c 



G 1 Ci C 1 C 2 

2 *— '1 »— '2 *--^2 



and recalling the fact that C is full-rank, we have that Cf C 2 -Cf Ci(Cf Ci) -1 Ci = Cf P Cl C 2 
is full rank. Similarly, one can show that (CfP Ca Ci) exists. Also, we have |C ff C| = 

CHr>± C* I \(-*H(-* I \f~fH-p± /-< I \(~^Hr^ I 



Define C = [P C2 Ci C 2 ], so that 



u 



C"C=[P£ 2 d C 2 ] [P C2 C X c 2 ] 



CfP C2 Ci 

cfc 2 



Consequently, we can write 



c 2 



C 2 (C 2 C 2 J C 2 — c 







Q (CfC.)- 1 



C 



fl 



(67) 



5 Note that ffS^ G^GSs 3> B is actually diagonal, and hence also block diagonal. 
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and 



c 2 



I-Pc 2 



C(C H C) C H -P C2 



-i 



(CfP^C,)- 1 Q 




lli 



(68) 



In the following we show the achievablity of R\ in §19) . The achievablity of _R 2 is obtained 
in a similar manner. Since V\ and V^ in Theorem [3] are independent, from (14) we have 



R x = I{V i; H(V! + V 2 ) + ZO - I(Vi; GFi + Z 2 ) 



logll + HSH^j -log 
log 1 1 + SH^Hl - log 



I + H(S2P C2 S2)H 



1-1 



log 



I + G(S2P^ 2 S2)G 



H 



I + S5P C ,S5H W H 



Hi 



-log 



I + S5P^ S^G^G 



Recalling (53), we have 



|I + SH H H| = l(C^C)- 1 ! • |Ai|-|A 2 | 

= |(CfPi a Ci)- 1 |.|(CfC 2 )- 1 |.|Ai|.|A 2 | , 



where we used Remark [5] to obtain ( f70] ). From ( |67| ), we have 



log 



I + S2P C2 S2H"H 



I + C 



1 + 



1 + 



1 + 





2 °2 

i) 

2 *^2 





(C^C^" 1 



(c^c,^ 1 



iH 



-H 



Hn-H 



C"C 



(c^c,)- 1 



Q Q 

Q (C^C,)- 1 



Ai 

A 2 

A! 
A 2 

I N^ Ai 

1 A 2 

A! + N^A 2 N N H \ 2 
A 2 N A 2 



C 



-1 



C^C-C^C 



I Q 
N I 



C^C 



-C*C 



(69) 



(70) 



(71) 
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(CfC 2 )- 1 A 2 N (CfC 2 )- 1 A 2 
CfCol-^ol = kef Car 1 ! • IA. 



j 2 V^2 



2 ^2) 



(72) 



where in (71 ), N = (Cf P^ 1 C 2 ) * Cf P^P^Ci, and we used the fact that 



C^C 



Similarly, we have 



(CfPjkCO^CfPc 



(C 2 Pci^i 



"CfP^ 



[Pc 2 Ci C 2 ] 



I Q 
N I 



log 



I + S5P^ 2 S5G^G 












1 + 


"(CfP^d)" 1 o" 




" I N H ~ 




I Q 




Q Q 




Q I 




N I 


1 + 


"(cfp&d)- 1 o" 




I + N^N N H 




Q Q 




N I 


I" (CfP^CO^CI + N^N) (CfP^d)" 1 ^ ] 











I 







c w c 



c w c 



KCfP^CO^I-ll + N^Nl 



(73) 



Subsituting (70), (72) and (73) in (69), we have R x = max(0, log |Ai| -log ll + N^Nl), which 



completes the proof. 

Appendix E 
Proof of Lemma[2] 

We want to show that WH fl HW and WG fl GW commute, where W = (U H U + G H G)-^ 
Let the invertible matrix C and diagonal matrix A respectively represent the generalized eigen- 
vectors and eigenvalues of (WH^HW + I ; WG ff GW + I), so that 



C H [WH fi HW + I] C = A 



(74) 
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C H [WG^GW + 1] C = I 



(75) 



Adding (74) and (75), we have 



C H [W(H H H + G H G) W + 2 1] C = 3 C^C = (A + I) 



from which it results that C must be of the form |[l2l 



C = -^* W (A + I)a 



(76) 



where $ w is an unknown unitary matrix. In the following, as we continue the proof, <fr w is 
characterized too. 



Substituting (76) in (74) and (75), it is revealed that the unitary matrix <fr w represents the 
common set of eigenvectors for the matrices WH fl HW + 1 and WG fl GW + 1, and thus both 
matrices commute. In particular, 

$^ [WH^HW + I] $ w = 3 A(A + I)" 1 = 3 (A" 1 + I)" 1 
^ [WG H GW + I] $ w = 3 (A + I)" 1 . 
Consequently, S x and S 2 are diagonal: 

$£ WH^HW $ w = 3 (A" 1 + I)" 1 - I = Si 
$£ WG fl GW $ w = 3 (A + I)" 1 - I = S 2 . 



(77) 



It is interesting to note that, since $^WH a HW$ w >z Q and $;WG a GW$ w >z Q, we 
have §1 r< A < 21. 



Appendix F 
Proof of Lemma[3] 

We first consider the generalized eigenvalue decomposition for 

(T^H"HT W + I , T^G H GT W + I) 



(78) 



April 18, 2013 



DRAFT 



30 



where T w is given by ([24]) and 



Cl [T^H H HT W + I] C w = A w 

C;[T^G H GT W + I]C W = I. 



Using (21 ), and noting that <fr w is unitary and P is diagonal, a straightforward calculation yields 



7 n 



C w [SxP + I] C w = A, 



7 H 



C w fS 2 P + IlC, 



I 



where £i and £ 2 are respectively (diagonal) eigenvalue matrices ofWH H HW and WG ff GW, 



as given by (21 ). Thus, C w is diagonal and is given by 



:s 2 p + i) 



(79) 



Consequently, we have A w = (£ 2 P + I) l (EiP + I). 

Let an, o 2 i and pi represent the i th diagonal elements of Si, S 2 and P, respectively. We 
note that for any pi, (1 + oiiPi)/(l + 02iPi) > 1 iff &u > cr 2 i. Thus, based on the argument that 
we made after Lemma [2] the first p diagonal elements of A w represent generalized eigenvalues 
greater than 1. Letting 



Pi Q 
o P 2 

where Pi is p x p and P 2 is (n t — p) x (n t — p), we have: 



(80) 



A, 



^Pi+ir'osipPi+i) 







1-1 



E^Pa + ir 1 S lp -P 2 + I 



lpJ 



Aiw 
Q A 2w 



(81) 



where T, ip and Ti ip (i = 1,2) are given by (22). Consequently, (79) can be rewritten as 

C w = [Ciw C 2w J = 



S 2p Pi+I) 











S 2/5 P 2 + 1) 



(82) 



From the argument before Lemma [3] for any diagonal 
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P >z 0, linear precoding is an optimal solution for the BC under the matrix power constraint 
S w = T W T^, where T w is given by (24). More precisely, from Theorem [Tj X — V\ + V 2 is 
optimal, where V\ and V 2 are independent Gaussian precoders, respectively corresponding to W\ 
and W 2 with zero means and covariance matrices K* w and S w — K^ w , where K^ w is given by 



K 



tw 



Ow^\ 



(Cl w Ci w ) 










CHq 2 
w°w 



(83) 



and C w is the generalized eigenvector matrix for 



(sIh^ 



rio w ~r A 1 Sw" Cjrovv t i- 



lit 



(84) 



We note that there exists a unitary matrix * for which S w = T W * H (H), where S w = T W T^. 
We also note that, from Remark |4J C w = *C W and C^C W = cf C w . Thus, K^ w can be 
rewritten as 



TV-* _ rp p 
xv (w J -w v -'w 



*H 



(C^C^)- 1 




7^ H r r H _ rp 
^w X w X w 



I Q 




1.H 



W$ W P2 



I Q 





w 



W$, 



Pi Q 








(85) 



(86) 



where (85) comes from (82), and (86) comes from (80). Consequently, S w — K* w can be written 

as 



K 



tw 



Trr\H 
w J- w 



K 



tw 



W$ W P$^W - K* w 



W$, 



Q Q 
P 2 



*™W 



(87) 



From ( [86] ) and ( [87] ), under the matrix power constraint S w given by ( [23] ), the optimal linear 
precoding is X = V\ + V 2 , where precoding signals V\ and V 2 are independent Gaussian vectors 



the optimal precoder can be represented as X = W<E\ 



with zero means and covariance matrices given by ( [86] ) and ( [87) , respectively. Alternatively, 



vi 



, where precoding signals V{ 
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and V2 are independent Gaussian vectors with zero means and diagonal covariance matrices 
respectively given by Pi and P 2 . In both cases E{XX H } = S w , and the same secrecy rate 
region is achieved. 
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